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Abstract 



The temperature dependence of the kinetic coefficients is obtained in the non- 
perturbative region with the help of Green-Kubo-type formulae in the model of 
massive gluon gas motivated by numerical results from simulations of lattice 
QCD. The entropy production rate is estimated using scaling hydrodynamics. It 
is shown that the increase in the viscosity coefficients leads to entropy generation 
in heavy ion collision processes which could be big, especially for temperatures 
close to the critical one. 



^Department of Physics, University of Cape Town, Rondebosch 7700, South Africa 
^Department of Physics, Saratov State University, 410071 Saratov, Russia 

■^Bogolubov Institut for Theoretical Physics, Academy of Science of Ukraine, 252143 Kiev, Ukraine. 



Discussions of the forthcoming projects for ultrarelativistic heavy ion coUisions 
at RHIC and LHC require, with necessity, comprehensive estimates of the space-time 
development in those reactions. Nowadays, it is commonly believed that central heavy 
ion collisions pass through several stages. Describing these in terms of equilibrium 
processes when all memory of prior history has been destroyed, it is natural to consider 
the entropy to set up the equilibration time scale and to separate a preequilibrium stage. 
In fact, that stage is producing not only the maximum entropy attained [1,2], but a lot 
of uncertainties in the initial and boundary conditions for the subsequent hydrodynamic 
expansion. It has usually been taken as isentropic, however it became clear that this 
stage may be complicated by the dissipative processes generating entropy, together 
with a possible phase transition (or transitions) when it happens to be of first order 
[3,4,5]. The aim of the present study is therefore to explore the entropy generation 
just at this stage although we understand that the freczc-out stage, where the system 
is made up of frcc-strcaming final particles, could also add to the entropy [6,7]. 

In order to estimate properly the dissipative effects as well as the dynamics [8] of the 
QCD phase transition we need to know the behaviour of the kinetic coefficients (KC) 
over a wide range of temperatures, including the ones close to the phase transition point 
where non-perturbative effects are dominant. The calculations done previously suffered 
from the unjustified extrapolation of the asymptotic behaviour, found pcrturbatively, 
to the critical region [9,10,11]. The recent estimate of the shear viscosity coefficient, 
using a model for the contribution of the nonperturbative region, has demonstrated 
that the behaviour in the critical region is very different from the standard one and 
the amount of entropy generated in the region close to this temperature appears to be 
substantial [12,13]. 

In those calculations we have exploited the so-called momentum "cut-off model" 
motivated by a special analysis of the numerical results of lattice Yang- Mills field 
thermodynamics [14,15]. In spite of the fact of providing a good fit to the data, this 
model is very much phenomcnological and doesn't interpret why the low-momentum 
modes are removed nor does it answer questions in terms of conventional conceptions 
of spin systems which are very indicative and conclusive at least for the Monte Carlo 
analysis of lattice pure gluodynamics. Here we are dealing with the model of an ideal 
gas of excitations where the effect of interactions in the plasma is provided by the 
temperature dependence of their effective mass [16]. In a sense it is rather similar to 
the "cut-off model" as now the medium (plasma) properties suppress low momentum 
excitations too, since for M{T) increasing cxp(— + p'^ /T) cxp{—M/T) at small 
momenta p — > 0. But it has more advantages as was argued in a recent analysis [17], 
the most important of them is that it leads to a perfectly detailed description of high 
precision SU(2) pure gluodynamics lattice data just in the weak coupling limit. 

However, as to the KC calculations, this model is more involved as in what follows 
we need to deal with massive scalar theory. By the way, this latter approximation 
is also in line with the model of "massive" gluons where these waves are nevertheless 
considered to be only transversal [18]. Combining our calculations of KC's in A0^- 
theory with the temperature dependence of the gluon mass extracted from lattice 
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Monte Carlo data, in particular for the mass gap in pure gluodynamics [19], we are 
able to estimate the entropy generated based on linear hydrodynamics and to explore 
the applicability of this approach to the evolution of gluon systems. 
As our basic hydrodynamical equation we take 

* + i±p_4 = o. (1) 

with Bjorken initial and boundary conditions [20]. We fix the equation of state (p = 
p(ro),£ = £(to) are the initial pressure and energy density respectively), and taking 
the initial conditions at a time Tq ~ 1 fm. The dissipative term in Eq.(l) contains 
the factor x — {^I'^'Hs + 'Hv) with ry^and r]y as transport coefficients of shear and bulk 
viscosities. We should as well take into account that t > otherwise we deal 

with the unrealistic picture of gluon gas contraction (see [5,11]). 
The total entropy of the system is defined as [5,10] 

j dai's'' = j dys{r)r , (2) 

where — su'^, s{t) = [£{t) +p(r)]/T is the local entropy density and y is the 
hydrodynamic rapidity (tanhy — x/t). Then Eqs.(l) and (2) give us a simple formula 
to estimate the entropy production in expanding gluon gas 

dS_ ^ r dTxjr) 
dy J tT{t) 

here %(''") = T{t) being the solution of Eq.(l). In order to solve Eq.(2) 

we need the temperature dependence of the KC's and the equation of in the whole 
temperature interval, including the nonperturbative critical region. 

For pressure and energy density calculations we use the familiar expressions [21] 

oo 

<T) = ^ I dkk'E{k)n{E{k)), (5) 



where E{k) = ^k^ + w?g is the relativistic energy; n{x) = [Gxp{(3x) — is Bose's 
distribution function; (3 is the inverse temperature and g is the degeneracy factor. In 
the ultrarelativistic case (T >> uig) it leads to the Stcphan-Boltzmann (SB) law (see 
Fig.(l)). However, if the temperature is close to Tc, the gluon mass rrig is increasing, 
(at least for SU(2)-gluodynamics) and becomes too large to be negligible and Bose's 
distribution function may be replaced by Boltzmann's {n{x) = exp(— /3a;)). Then eqs 
(6), (7) may be taken in the following form 

P{T) = ^-^Jdz[z'- aT' = |^ «^ Mo^), (6) 

a 
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a 

where Ki{a) are modified Bessel function and a = nigP. It is a well-known fact 
that the asymptotic expansion of the modified Bessel functions is given by K{a) ~ 



ex.p{—a)^7r/2a for large a, i.e. for T ~ Tc, this way both pressure and energy become 
small close to Tc. In the case when a — > we recover the Stephan-Boltzmann law. 
This temperature dependence can fit even the SU(3) Monte-Carlo lattice data with 
constant gluon mass [18] or including finite jump of mass around Tc. As to SU(2)- 
gluodynamics, for which a much more elaborated Monte Carlo analysis exists [17], we 
used the following parametrization 



m, = moTc (^) 



with rriQ = 1.83 and q = 0.4. An important difference just reflects our understanding in 
behaviours of first- and second order phase transitions as seen in lattice Monte-Carlo 
simulations [19]. To calculate the KC's of shear and bulk viscosities we use well- 
known relations obtained within a formalism based on Kubo - type formulae for the 
- thermofield theory [10,11] (The analogous expression for rjs was also obtained for 
the case of vector fields [22]) 

^. = ^^2,1, (8) 

Vv = ^{/2,i - Gc^hfi + 9c,%_i}, (9) 
where = dp/de is the square of sound velocity, and the integrals Im,n ^-re defined as 

W = 2 / n{p) [ 1 + n{p% (10) 

here r(p) is the damping rate of quasiparticle excitation (it is assumed that V(5 « 1). 
For the scalar theory in one- loop approximation is 

A^(27r)^ r ^ ^ ^ 
r(P) = OA ^ \ ^ / d%d%d%d{p + pi-p2-p3)il + ni)n2n3 (11) 

with the following designations d^p = [2{27r)^E{p)]~^d^p, in both Eqs.(ll) and (12), 

The gluon gas viscosities may be obtained from Eqs.(8)-(11) by the standard procedure 
of changing [11,23] 

— ^ c*327r2«2 In a;\ c* = 20 60, (12) 
where (for the value of c*, see also [24]) 



6% 



-1 



ll/2iVln(M7A') , (13) 
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and is a number of colours, = | < ^ g^j^^j < > is the thermodynamicaly 
averaged squared momentum of the gluon field [21]. The degeneracy factor g — [d — 
1){N'^ — 1) is absent in final result, for the numerator of Eq.(12) must contain it as well 
as the damping rate in the denominator. This phenomenological estimate can also be 
justified by the fact that in the lowest order of interaction the cross sections for gluons 
and scalar particles have a similar momentum dependence. This procedure (Eq.(13)) is 
fair only for small enough interaction constant. The model under consideration brings 
us to the following temperature dependence for the factor from Eq.(15) 

4 Jd^pn{E{p)) K^ia)- ^ > 

In the ultrarelativistic case where a — > we obtain the conventional result M ~ 4T. 
When T — >■ Tg, and a — > cxd we have ~ 4T^q;. This means that the coupling 
constant ag remains still small even when the temperature is close to T^. This lucky 
fact has been met already in the so-called cut-off model [14,15] that interprets Monte- 
Carlo data as well as our model does. It can be explained evidently by the fact that 
the application of both gluodynamic models takes the contribution of long wavelength 
excitations away. 

We now proceed to calculate the damping rate (Eq.(ll)). After some straightfor- 
ward integrations using the 5 - function we have 



192E(p)n{E(p)) J E{p,) 
where 

h = 2np-' I dy exp(-/?n) e{zo{y) - 1) e{zo{y) + 



My) - ^!^^:±S, (16) 

here we need to solve an inequality 

- 1 < z^iy) < 1, (17) 

where 



Q = y^H^ + v^p2 + ml (18) 

and K — ^Jp^ +Pi + 2ppi cos(^), 6 being the angle between p and pi. It then leads 
us to the expression 

(19) 
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with the values of Q and K satisfying the inequality 
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After analysis of Eqs. (19) and (20) we obtain for Eq.(15) 



^2-/^2 > ^20) 



These relations are complicated to integrate, in order to get an approximate analyt- 
ical form we can calculate the angle integral from Eq. (21) in the p = case. This gives 
a minimal value for the damping rate and a maximum one for the kinetic coefficients. 
Here we also take into account that the main contribution in Eqs. (8)- (10) is connected 
with long wavelenght excitations. Thus, in this approximation the damping rate looks 
as ^ 

r(p) = I dz{z- a) (e- + e-^^) . (22) 

a 

It is evident that the momentum dependence is very simple in this equation. This 
allows us to represent the integral from Eq.(lO) in the following form 

rp2{n-m)-2 

^ 2 7r2f(T) "^""'^ ^^^^ 

with 

CXD 

Jm,n = / dzz^-^"" {z^ - a^f-^[exp{-z) + exp(-2^)] (24) 

a 

After integration we obtain, for some specific cases 

J2,i = ^ [(2a)^K:,(a) + a^K^(2a)], (25) 

Ji,o = 3a^ [- Ks{2a) + 5 ^3(0;) + 7 i^2(2a) + aK2{a)l (26) 

o 4 

Jo,-i = a' [60X2(0;) + ^ K2{2a) + 27 a Ki{a) + ^ Ki{2a)+ 

4 8 

+ 6a'^ Ko{a) + -a^Ko{2a) + a^Ki{a) + — Ki{2a). (27) 
The result for the damping increment can be written as 

f(r) = E{p)T'r{p) 

\ 2 —3 

= ^|^«2(T)[2ir2(a)-2i^i(a) + K2(2a)-i^i(2a)] (28) 
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For the calculation of the bulk viscosity we also need to know the expression for 
the velocity of sound. With the help of Eqs. (6) and (7) we obtain 

' 3 + A 

with 

Here we see that, for a small, A ~ a, and the sound velocity vanishes. In the a — > 
case Eq.(28) gives the common value of = 1/3, A — > (see Fig. (2)). 

Eqs. (24-28) together with Eqs. (8, 9, 10, 12-14) determine the temperature depen- 
dence of the KC of the gluon gas. The asymptotic behaviour a — > gives results 
analogous to Ref. [11] that is somewhat large, for, here we have a maximal estimate 
of KC (see below Eq.(21)). The temperature dependences of the KC are presented in 
Fig. (3). It is evident that the KC increase considerably in the temperature region close 
to Tc- This leads to a big deviation of the solution of Eq. (1) from the scaling one, and 
to the so-called critical delay of the evolution of the system near Tc. 

The rate of entropy production as a function of proper time is depicted in Fig. (5) . It 

is evident that the entropy of the system increases rapidly. The analogous calculations 
for entropy production rate for T^- dependency of KC brings to approximately 20 per 
cent increasing of entropy in gluon gas cooling process. 

The authors are aware of the fact that the obtained results are model dependent. 
The increase of the kinetic coefficients makes the application of linear hydrodynamics 
for the description of quark gluon gas questionable. Evidently, taking into account other 
dissipative mechanisms will lead to finite values for the KC and diminish somehow the 
entropy production. However the results obtained indicate that a realistic picture of 
the evolution of the system under consideration can differ a lot from the scaling one, 
especially in the phase transition region and demands thorough examination. 
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Figure Captions 



Figure 1. Dependence of the scaled energy density and pressure on the temperature. 

Figure 2. Dependence of the speed of sound on the temperature. 

Figure 3. Evolution of the temperature with proper time. 

Figure 4. Dependence of the shear and bulk viscosities on the temperature. 

Figure 5. Entropy generation as a function of proper time. 
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